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ABSTRACT 


SEIMR/R-S corresponds to a generalized mathematical model of pandemics that 
enhances traditional, aggregated simulation models when considering inter¬ 
regional impacts in a macro region (conurbed); SEIMR/R-S also considers the 
impact of modeling the population divided into socio-demographic segments based 
on age and economic stratum (it is possible to include other dimensions, for 
example: ethnics, sex,...). 

SEIMR/R-S is the core of the SEIMR/R-S/OPT epidemic management optimization 
model that determines optimal policies (mitigation and confinement) considering 
the spatial distribution of the population, segmented socio-demographically. The 
formulation of SEIMR/R-S/OPT is presented in another "paper" describing its 
implementation in GAMS and AM PL, and the implementation for the case of 
Bogota. 

SEIMR/R-S can be understood and used by any epidemiologist, and/or physician, 
working with SIR, SEIR or similar simulation models, and by professionals working 
on the issue of public policies for epidemic control. 


1. EPIDEMIC & CONTROL POLICIES MODEL 

The SEIMR/R-S is a detailed epidemic model that is the result of integrating the SIR, SEIR and 
SEI3RD model; in these standard models the population is grouped in only one homogenous group. 
SEIMR/R-S extends the modeling to a multi-segment-sociodemographic multi-region system. 

SEIMR/R-S model describes the epidemic with following states: 

S Susceptible: initially covers all population that potentially can be infected (SU) 

E Exposed: Population that has been infected and are in an incubation (latency) period (EX). The 

model SIR does not include this state. 

IM Multi-Infected: Population that has been infected and has active the pathogen in different states 
of development (10,11,12,..., IN). The active infected states are ordered according to the severity 
of the infection. The modeled SIR and SEIR consider only one infected state. For convenience, the 
last state is called "IN" 

R Recovered: Recovering population (RE) 

R-S is related with the Region-Segment model that considers multiples regions where live people 
classified in multiples socio-demographics segments. 

1.1. CONCEPTUAL FRAMEWORK 

The following documents have been referenced and used for writing the following numerals: 

■ Brauer, F. and Castillo-Chavez, C. (2000). Mathematical Models in Population Biology and 
Epidemiology. Springer, New York, 2000. 

■ Cai, Y., Kang, Y., Banetjee, M., & Wang, W. (2015). A stochastic SIRS epidemic model with infectious 
force under intervention strategies. Journal of Differential Equations, 259(12), 7463-7502. 

■ Carcione, J., Santos J. E., Bagaini, C. and Jing, Ba. (2020) A Simulation of a COVID-19 Epidemic 
Based on a Deterministic SEIR Model. 

■ Erdem, M., Safan, M., & Castillo-Chavez, C. (2017). Mathematical analysis of an SIQR influenza model 
with imperfect quarantine. Bulletin of mathematical biology, 79(7), 1612-1636. 

■ Eubank, S., Eckstrand, I., Lewis, B., Venkatramanan, S., Marathe, M., & Barrett, C. L. (2020). 
Commentary on Ferguson, et al.,"Impact of Non-pharmaceutical Interventions (NPIs) to Reduce 
COVID-19 Mortality and Healthcare Demand". Bulletin of Mathematical Biology, 82,1-7. 
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■ Ferguson, N., Laydon, D., Nedjati Gilani, G., Imai, N., Ainslie, K., Baguelin, M.,... & Dighe, A. (2020). 
Report 9: Impact of non-pharmaceutical interventions (NPIs) to reduce COVID19 mortality and 
healthcare demand. 

■ Hethcote, H. W. (2000). "The Mathematics of Infectious Diseases." SIAM Review 42.4 (2000), pp. 
599-653. 

■ Huang G. (2016). Artificial infectious disease optimization: A SEIQR epidemic dynamic model-based 

function optimization algorithm. Swarm and evolutionary computation, 27, 31-67. 

https://doi.ora/10.1016/i.swevo.2015.09.007 

■ Jing, W., Jin, Z., & Zhang, J. (2018). An SIR pairwise epidemic model with infection age and 
demography. Journal of biological dynamics, 12(1), 486-508. 

■ Kermack, W.O. and Me Kendrick, A. G.,(1927) "Contributions to the Mathematical Theory of 
Epidemics, Part I". Proc. R. Soc. London, Ser. A, 115, 1927, 700-721. 

■ Liu M. and Liang, J. (2013). Dynamic optimization model for allocating medical resources in 
epidemic controlling. Journal of Industrial Engineering and Management JIEM, 2013 - 6(l):73-88 - 
Online ISSN: 2013-0953 

■ Massed, E., Bruattini, N.M., Coutinho, B. A. F. and Lopez, F. L., (2007). "The 1918 Influenza An 
Epidemic in the City" of Sao Paulo", Brazil, Medical Hypotheses, 68, 2007, 442-445 

■ Mejia Becerra, J. D. et. al. (2020). "Modelacion Matematica de la Propagacion del SARS-CoV-2 en la 
Ciudad de Bogota. Documento de Circulacion Informal 

■ Pang, W. (2020). Public Health Policy: COVID-19 Epidemic and SEIR Model with Asymptomatic Viral 
Carriers. Department of Mathematics and Statistics, McMaster University arXiv:2004.06311vl [q- 
bio.PE] 14 Apr 2020. 

■ Radulescu, A., & Cavanagh, K. (2020). Management strategies in a SEIR model of COVID 19 
community spread. arXiv preprint arXiv:2003.11150. 

■ Samsuzzoha, Md. (2012) . "A Study on Numerical Solutions of Epidemic Models". Doctoral Thesis, 
Swinburne University of Technology, Australia, 2012. 


The goal of epidemic control strategies is to reduce Ro. This can be achieved by reducing susceptibility or 
contact rates in the population or the infectiousness of infected populations. The potential effectiveness 
of medical intervention by varying the infectiousness of infected populations and nonmedical interventions 
by reducing the contact rates in the population have been examined. In medical intervention, use of 
vaccines and/or antiviral agents for case of treatment can increase the recovery rate and reduce the 
death rate. On the other hand, in nonmedical interventions, reducing population contact rates through 
social distancing and travel restrictions can reduce the impact on the transmission process. 




GOAL OF EPIDEMIC MODELING 


FIGURE 1. Goals of community mitigation for pandemic influenza 


mu 



Source; Adapted from: CDC. Interim pre-pandemic planning guidance; community strategy for pandemic influertza mitigation in the United States—early, targeted, 
layered use of nonpharmaceutical interventions. Atlanta, GA: US Department of Health and Human Services, CDC; 2007. httpsV/stacks.cdc.gov/view/cdc/l 1425. 


Control of an outbreak relies partly on identification of the disease parameters that lead to a significant 
reduction of the basic reproduction number Ro that may be function of several parameters of which y, 
the recovery rate for clinically ill and p, the transmission coefficient, are the most sensitive parameters. 
These two parameters can be controlled by medical intervention and nonmedical interventions. 
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The modeling of epidemics in a solidly developed area of scientific knowledge, widely studied based on 
simulation models. The master table of epidemic models shows some of the best-known models 


TABLE: MAE EMO 


STATE 

DESCRIPTION 

Reference 

Imple 

(COD EMO) 

(DIN EMO) 

(COM EMO) 

mented 

SIR 

Susceptibility (S), Infection (I) and Recovery (R) 

Kermack & Me Kendrick (1927) 
Jing (2018) 

YES 

SEIR 

Susceptibility (S), Exposure (E), Infection (I) and Recovery 
(R) 

Hethcote (2000) 

YES 

SEIRA 

Susceptibility (S), Exposure (E), Infection (I) and Recovery 
(R) 


NOT 

SEI3RD 

Susceptibility (S), Exposure (E), 3-1-1 Infection States (13), 
Recovery (R) and Death (D) 

Mejia Becerra et. al. (2020) 

YES 

SEIQR 

Susceptibility (S), Exposure (E), Infection (I), Quarantine (Q) 
and Recovery (R) 

Huang (2016) 

NOT 

SIRS 

Susceptibility (S), Infection (I), Recovery (R) and 

Susceptibility (S) 

Cai (2015) 

NOT 




NOT 

SEIMR/R-S 



YES 


This Part presents the mathematical formulation of the SIR and SEIR models. The SEI3RD is presented 

in the Part IV; Epidemic Optimization Modei Appiied to Bogota City. 

Then, an epidemiological model is defined based on differential equations that explain the evolution of 
the process without human intervention. These differential equations can be established based on the 
population (number of people) who are in a certain "epidemic" state or based on the fraction of the 
population that is in that state. The epidemic models are nonlinear systems of ordinary differential 
equations, traditionally this equations system is solved using simulation models based in a discrete 
approximation for continuous derivatives, be it over time or space. There are many possible schemes. 
These models are used to analyze several widely discussed (predefined) scenarios and provide evidence 
on their effectiveness and are not oriented to get the optimal solution of a mix of control policies. 

The added value by mathematical programming approach is to convert simulation models into 
optimization models to be able to combine them with other mathematical programming models, following 
the principles of structured mathematical modeling that allows join multiple problems of mathematical 
programming in a single model. Based on the above, the formulation of the models is done by means of 
algebraic equations that represent how the epidemiological process evolves during the planning horizon. 

For the optimization epidemic modeling, the approach is based on multiple state chains that can be 
associated with semi-Markov chains; initially, it was proposed to model based on the approach of semi- 
Markovian processes (changing transition matrices over time) but such an approach brings multiple 
complications in the math formulation of probabilities. 
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After analyzing the implementation of main (most known) epidemiological models (SIR, SEIR), it was 
decided to directly model discrete versions of differential equations as they maintain direct connection 
with biological parameters, which facilitates the connection of these parameters with socio-demographic 
segments. 

Therefore, all epidemiological models considered should be formulated in one of the following terms. 

■ The time unit of the differential equations is one day. 

■ The states contain the fraction of the population in each state. 

■ The time of the optimization model may be divided in periods of multiple days (one week, seven 
days). In this case, the integration of the differential equations must be made using calculated 
parameters. 

The epidemic states are showed in the master table MAE_STA that describes the states. The models 
will be implemented using this nomenclature. The table includes the symbol used in the original models 
and the code using in the information system to reference the state. 


EPIDEMIC STATES - TABLE: MAE STA 

MODEL 

SYMBOL 

EPIDEMIC 
STATE 
(COD STA) 

DESCRIPTION 

(DIN_STA) 

COMMENTS 

(DLIN_STA) 

S 

su 

Susceptible 

Population 

Those individuals who have not been exposed to the pathogen and are 
susceptible to being infected by it. 

E 

EX 

Exposecj Population 

Those individuals who are in the latency state; that is, they have been 
inoculated by the pathogen but are not yet infectious 

I 

IN 

Infectecj Population 

In SIR and SEIR models is infected population. It must be the most 
critical state for infected people; this is important for models that have 
more than one epidemic states to describe the infection process. 


lU 

Infectecj Population 
Unatten(je(j 

Individuals that require hospital attention and do not received hospital 
ty attention 


ED 

Epidemic Dead 

Individuals who fail the infection and die. 


CD 

Collateral Dead 

Individuals who die by during the epidemic by reasons different to the 
epidemic. Minimize this population is the goal of the optimization when 
the epidemic in management isolated, ignoring the economic impact. 


ND 

Natural Dead 

Individuals who die by other reason different to the epidemic 


NP 

New Population 


R 

RE 

Recovered 

Population 

Those individuals recover from infection, having developed antibodies. 

In most of the models they cannot be re-infected. 


EPIDEMIC STATES - TABLE: MAE STA 

MODEL 

SYMBOL 

EPIDEMIC 
STATE 
(COD STA) 

DESCRIPTION 

(DIN_STA) 

COMMENTS 

(DLIN_STA) 

STANDARD STATES MBC-SEI3RD MODEL 


ED 

Epidemic Dead 

Individuals who fail the infection and die. 

E 

EX 

Exposed Population 

Those individuals who are in the latency state; that is, they have been 
inoculated by the pathogen but are not yet infectious 

lo 

10 

Asymptomatic 

Infectious 

Those individuals in the population who have been inoculated by the 
virus are infectious but have not developed symptoms. Those infected 
in this state rarely learn that they have been infected. 

Ii 

11 

Moderate Symptoms 
Infectious 

Those individuals in the population who are infectious and have mild or 
moderate symptoms. They are those who can be given management of 
the disease at home. 

h 

12 

Severe Symptoms 
Infectious 

Those individuals in the population who are infectious and have severe 
but not critical symptoms. Individuals present in this state require 
hospitalization. 

h 

IN 

Critical Symptoms 
Infectious 

It must be the most critical state for infected people; this is important 
for models that have more than one epidemic states to describe the 
infection process. In SIR and SEIR models is infected population 


ND 

Natural Dead 

Individuals who die by other reason different to the epidemic 


NP 

New Population 


R 

RE 

Recovered 

Population 

Those individuals recover from infection, having developed antibodies. 

In most of the models they cannot be re-infected. 

S 

SU 

Susceptible 

Population 

Those individuals who have not been exposed to the pathogen and are 
susceptible to being infected by it. 
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EPIDEMIC STATES - TABLE: MAE STA 

MODEL 

EPIDEMIC 
STATE 
(COD STA) 

DESCRIPTION 

COMMENTS 

SYMBOL 

(DIN_STA) 

(DLIN_STA) 

ADDITIONAL CAPACITY STATES MBC-SEI3RD MODEL 


IU 2 IU 3 

Infected Population 

Individuals that require hospital attention and do not received hospital 


Unattended 

ty attention in state I 2 , 13 , ... 


CD 

Collateral Dead 

Individuals who die by during the epidemic by reasons different to the 
epidemic. 


The measurements used must be equal for all models 


Measurement Unit Master Table: MAE UND 

Measurement Unit 

Description 

COD UND 

DES UND 

1/peo-day 

1/ persons-day 

fpo/day 

Fraction of population per day 

peo-day 

Persons-day 


One of the main limitations of the traditional approach is to assume that the entire population is 
homogeneous with respect to its epidemiological behavior. It is well known that the epidemic manifests 
differently in each socio-demographic stratum and that the composition of socio-demographic segments 
depends on each region. 

In order to enhance the model to be useful in real cases, it is assumed that there is a different pandemic 
(because it has different parameters) for each pair <region, demographic-segment These hypotheses 
may vary according to each case study. In this case, reference has been made to the data used to control 
the epidemic used in Bogota. It should be noted that the parameters of each epidemiological model vary 
in quantity but cannot vary in the form of calculation, since they are parameters as the same case, which 
is studied with different mathematical models. 

The models are studied under the hypothesis of a homogeneous population in a region, then the epidemic 
is assumed to be particular to each duple <rg,ss> and the equations are formulated depending on 
<rg,ss>. The advantage of this approach will be visualized when the epidemic model is coupled with the 
management of health resources and control policies, which can be individualized for each duple <rg,ss>. 

The model parameters can be grouped by the original source of variation, these sources are: 

■ Pathogen: characteristics of the epidemic due to the pathogen 

■ Age: It is typical for recovery/worsening times (rates) and probability of recovery to be a function of 
age. 

■ Economic stratum: influences the epidemic by means of the intensity of contact, product of the 
number of contacts, the duration of contacts and the closeness, these variables may also be a 
combined function of age and economic. 

Additionally, may be considered people coming for the exogenous systems (out of the microregion) to 
the region. 

The indexes used in the modeling are presented in the next table. 


INDEXES 

Index 

HEA 

Short Description 
(Entities) 

Long Description 

ag 

Age 

Age 

mr 

Macro-region 

Macro-Region 

rg, ro, rd 

Region 

Region (Basic Territory Unit) 

ss 

Social Segment 

Socio-Demographic Segment 

St, si 

Epidemic State 

Epidemic State 

t,q 

Period 



In the Bogota case, the socio-demographic segments are a combination of age with an economic stratum. 
The biological parameters depend on age. 
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2. GENERAL SIMULATION EPIDEMIC MODEL 

Below is presented an aggregate model of epidemic that is the result of integrating the SIR model, SEIR 

and SEI3RD; in these standard models the population is grouped in only one homogenous group. 

SEIMR/R-S extends the modeling to a multi-segment-sociodemographic multi-region system. 

SEIMR/R-S model describes the epidemic with following states: 

S Susceptible: initially covers all population that potentially can be infected (SU) 

E Exposed; Population that has been infected and are in an incubation (latency) period (EX). The 

model SIR does not include this state. 

IM Multi-Infected: Population that has been infected and has active the pathogen in different states 
of development (10,11,12,..., IN). The active infected states are ordered according to the severity 
of the infection. The modeled SIR and SEIR consider only one infected state. For convenience, the 
last state is called "IN" 

R Recovered; Recovering population (RE) 

R-S is related with the Region-Segment model that considers multiples regions where live people 

classified in multiples socio-demographics segments. 

The next table shows the relation between models and epidemic states. 


Model 

SEIMR/R-S Model Epidemic States 

Standard 

Extended 

Capacity 

SU 

EX 

lo 

II 

l2 

... 

In 

RE 

NP 

ED 

ND 

lU 

CD 

SIR 

X 






X 

X 






SEIR 

X 

X 





X 

X 






SEI3RD 

X 

X 

X 

X 

X 

X 

X 

X 


X 




SEIMR/R-S 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


The biological parameters used in the are described below (it includes the "standard" original formulation 
and the OPTEX formulation). 


BIOLOGICAL PARAMETERS (READ) 

Parameter 

OPTEX 

Parameter 

Table 

Description 

Measure 

Unit 

Model 

PATHOGEN BIOLOGICAL PARAMETERS - SIR/SEIR 

5 

DELT 

PBIO GEN 

Contact Intensity - Exogenous Parameter 

peo-day 

SIR 

CO 

PTRA 

PBIO GEN 

Probability of transmission per contact intensity (infectivity) 



Y 

GAGA 

PBIO GEN 

Recovery rate for clinically ill (independent of age) 

fpo/day 

SIR 

Yst 

GAGA 

PBIO_GEN 

Recovery rate for clinically ill (independent of age). It depends 
of St for SEIR3D 

fpo/day 

SIR 


ZETA 

PBIO GEN 

Epidemic death rate (independent of age) 

fpo/day 

SIR 


MIUN 

PBIO GEN 

Natural mortality rate (independent of age) 

fpo/day 


K 

KAPP 

PBIO GEN 

The latency period of the virus before developing 

day 

SEIR 

AGE BIOLOGICAL PARAMETERS - REGIONAL - SEGMENTED 

|iag 

MlUUag 


Epidemic mortality rate (complication rate, age dependent). 

fpo/day 

R-S 

yaag 

GAMUag 


Recovery rate for clinically ill (age dependent) 

fpo/day 

R-S 

EXOGENOUS SYSTEM PARAMETERS - SEIR 

Xe 

LAMErg,ss 

UBT SDS 

Exposed rate coming from the exogenous system 

fpo/day 



LAMSrg,ss 

UBT SDS 

Susceptible rate coming from the exogenous system 

fpo/day 



LAMIrg,ss 

UBT SDS 

Infectious rate coming from the exogenous system 

fpo/day 


Xr 

LAMRrg,ss 

UBT SDS 

Recovered rate coming from the exogenous system 

fpo/day 


BIOLOGICAL PARAMETERS - SEIR 


BEQUrass 

1 1 

Parameter p for confinement control policies (SIR) 




The calculated biological parameters used are 
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BIOLOGICAL PARAMETERS (CALCULATED) 

Parameter 

Equation 

OPTEX 

Parameter 

OPTEX 

Equation 

Description 

Measure 

Unit 


1/k 

FHII 

1/KAPP 

Inverse virus latency period 

1 /day 

p 

55 X CO 

BETA 

DEIl X PTRA 

Inverse contact intensity x infectivity 

1 /fpo-day 


5-^ 

DEIl 

1/DELT 

Inverse contact intensity 


ILlUag 

SageAGS(ss) ILlag 

MISSss 

ZageAGS(ss) MlUUag 



YYag 

2ageAGS(ss) ILlag 

GASSss 

2^ageAGS(ss) GAMUag 



YILlss 

21ageAGS(ss) 
(Yag + ^lag) 

GAZEss 

GASSss + MISSss 




The calculus of parameters is analyzed in a posterior section, in the presentation of Bogota case. 

The general assumptions for standard epidemic models are: 

■ No vaccine exists 

■ The susceptible population is reduced through infection (moving to infective state). 

■ People who recovered after catch the virus will be insusceptible of it 

■ The population of infective class is increased by a fraction of susceptible individuals becoming 
infective. 

■ All other people are susceptible 

■ The population is homogenous 

■ The population of "critical" infective individuals is reduced by recovery from the disease. 

2.1. SIR; EPIDEMIC MODEL 

The SIR model is a basic model in epidemic modeling (Kermack and Me Kendrick, 1927). SIR process, 
starting with a susceptible host who becomes infected, the class of infection grow for the infected 
individuals to be able to transmit the infection to susceptible. When the infected individual is no longer 
able to transmit infection to susceptible individual, the infected individual is removed from the cycle of 
diseases transmission in the population. This model is based on the following assumption: 

Then, the basic SIR model describes the epidemic with three states: 

S Susceptible: initially covers all population that potentially can be infected (SU) 

I Infected: Population that has been infected (IN) 

R Recovered: Recovering population (RE) 

The diagram shows the behavior of S(t), I(t), and R(t) when they are normalized to total of population 
(TPOB) equal to 1. 

Susceptible-Infectious-Recovered (SIR) a I 

Epidemic Model -33AA/ 

- - - - - - DecIBIOMMrARE 



The basic SIR model is represented based on three differential equations based on proportions of people 
in each state (the ratio between the people in a state with the initial population TPOB). The 
measurements in blue. 

5S(t)/at (fpo/day) = - p (1/fpo-day) x I(t) (fpo) x S(t) (fpo) 
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5I(t)/5t (fpo/day) = p (1/fpo-day) x I(t) (fpo) x S(t) (fpo) - y (fpo/day) x I(t) (fpo) 

5R(t)/0t (fpo/day) = y (fpo/day) x I(t) (fpo) 

where S(t), I(t), R(t) represent the population of susceptible, infected and recovered individuals, 
respectively. Adding these equations, the following condition must be hold 


dS{t)/dt + dl{t)/dt + 5R(t)/5t = 0 

Additionally, SIR can be extended with other epidemic states for a more complete description of the 
system/epidemic: 

NP New population entering the system, as people from abroad who in many cases are the ones who 
cause the epidemic (NP). 

ED People who die due to the epidemic (these people die regardless of the management of the 
epidemic (D). 

ND People who die from natural death (N) 


ND and ED states should be included if it wants to account for the resources consumed by people who 
die, who are killed due the epidemic and due by causes other than the epidemic. The NP(t) state 
represent the proportion of people arriving from an exogenous system, may be births or people arriving 
from a foreign country/region. The value of NP(t) is a border condition with the foreign system over any 
value of t it is calculated taking as reference the initial population TPOB. This adjustment is made in all 
epidemic model included in the optimization model. 


STATE TRANSITION DIAGRAM FOR A MARKOV CHAIN MODEL 
SIR MODEL 


Wd 



STATE 

DESCRIPTION 

NP 

New Population 

SU 

Susceptive Population 

IN 

Infective Population 

RE 

Recovered Population 

ND 

Natural Dead 

ED 

Epidemic Dead 


I Smceptible | 

k 

I Infrcttoun I 
7 

I Rfco»crtd I 


dt 


=-psi 




dt'' 


ri 


Then, the differential equations must be adjusted: 

as(t)/at = - p X i(t) x s(t) + X NP(t) - x s(t) 

5I(t)/5t = p (I(t) X S(t) - (y + Q ] X I(t) + X NP(t) 

5R(t)/5t = y X I(t) - X R(t) + x NP(t) 

aD(t)/at = c X i(t) 
aN(t)/at = pN xS(t) + pN xR(t) 

where p'^ represents the natural mortality rate and A,®* the rates coming from the exogenous system to 
the state st. 
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If TPOB is the initial total population, and it is constant over time, NP(t)=0, the model meets the 
hypothesis that at all times 

S(t) + I(t) + R(t) + D(t) + N(t) = 1 
If NP(t) is different from zero the previous equation must be adjusted as 

S(t) + I(t) + R(t) + D(t) + N(t) = 1 + Jq^[o,t] aNP(q) Sq 

To simulate the process the border conditions at the beginning of the simulation horizon are: S(0), 1(0), 
R(0), D(0) and NP(t), for all t. 

Assuming NP(t) equal to zero, the ratio p = y/p is called the relative removal rate. Thus, dynamics of 
infectious depends on the following ratio: 


Ro = S(0) X y/p 

where Ro, called the basic reproduction ratio/number, is defined as the number of secondary infections 
produced by a single infectious individual during his/her entire infectious period. The role of the basic 
reproduction number is especially important. However, the following mathematical analysis describes how 
the basic reproduction number depends on the host population and the infected host. 

At time t = 0, dl/dt can be written as 

01/at = (Ro - 1) X y X 1(0) 

if Ro > 1 then ai/at > 0 and therefore the disease can spread; but if Ro < 1 then the disease dies out. 
Making mathematical manipulation it is possible to prove that the maximum number of infective at any 
time is 


TPOB (1 - p + p In [ p / S(0) ] ) 

It should be noted that the probability of transition becomes a dynamic variable that must be calculated 
by the mathematical model, for that reason the t index must be included. 

2.2. SEIR EPIDEMIC MODEL 

There are multiple versions of the SEIR model, below are some images summarizing some of the literature 
consulted. 


SOME SEIR MODELS 

Radulescu & Cavanagh (2020) 

^ = -pS(I+qE)/N 

^ = pS(I + qE)/N-SE 

f = 

dn 

Hethcote, H. W. (2000). 

+ vN + 6R 

at A 

dE ^TS , 

— = - (A* + cx + K)i. 

^ - (ac + a + 'y)l 

— kE ■+■ '■)I — fift — srt 

at 

aud 

S ^ E + I N 


Carcione et a. (2020) 

Pang, W. (2020). 
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SOME SEIR MODELS 

E = ^ 5 ^-(^ + £)E, 

/ = €£' — (7 + ^ -h q)/. 

jR = 7/ — fiR. 

^ = a'Sth + cc'^P^StEt - p'Et - 

dE, 

Liu and 

Liang (2013). 

dt 

— = p<k> S{t)I{t)-p<k>S{t-T)I(t- r) 
dt 

— = P<k>Sit-T)I(t-T)-id + S) lit) 
dt 

^=sm 

Vdt 

Grimm et al. (2020) 

d/f 


The classic SEIR model (Hethcote, 2000) describes the epidemic dynamics based on the transitions 
between four different compartments: susceptible (S), exposed (E), infectious (I), and recovered (R) 
individuals. The SEIR symbols are the same SIR symbols plus the parameter \|/ that represents the inverse 
of the virus latency/incubation period (k). 



The equations of SEIR model are the same as those of SIR model considering the following changes: 

■ Equation 5E(t)/5t is introduced 

5E(t)/5t = p X I(t) X S(t) - \|/ X E(t) + X^x NP(t) 

■ Equation 5I(t)/5t is modified 

ai(t)/at = ly X E(t) - (y + Q X I(t) + A,' X NP(t) 

2.3. SEI3RD EPIDEMIC MODEL 

SEI3RD was developed with the aim of simulating the transmission and evolution of acute infections. This 
simulation assumes that the pathogen causes an infection followed by lifelong immunity or death through 
multiple infected states. SEI3RD model extended the classic SEIR model to distinguish between several 
categories of infectiousness; for example: asymptomatic, symptomatic, moderate, severe cases, and ... . 
Being able to explicitly distinguish these different infected groups is important as they can greatly differ 
in terms of their underlying parameters as well as in terms of their behavioral response to public health 
interventions compared with the SEIR Model. 


11 




























m«TECNM 

mJ r^MEXlCO' 


The SEIMR/R-S General Epidemic. 

Theory Validation and Application in Mexican Regions 


DECISIONVIMRE 


STATE TRANSITION DIAGRAM - SEI3RD EPIDEMIC MODEL 



The transition between infected states follows this assumptions: 

■ A person can only be infected by one of the individuals belonging to one of the infected states (lo, 

II/ l2/ IS/ ■■■ / In^i 

■ Being inoculated by the pathogen, the individual passes to the group of exposed E. 

■ After the latency period, the person in the exposed state becomes asymptomatic infectious (lo state). 

■ Once the person is in 1st state, there are two possible outcomes: worsening clinical status (moving to 
mild/moderate/critical infected Ist+i) or recovery (R). 

■ If the person recovers, in any of the states of infection, they enter the absorbent recupered state (R). 

■ Similar logic applies for a larger number of states of infection. It is understood that if an individual of 
the last state of infection (13) worsens its clinical condition if it die (enters state E). 


BIOLOGICAL PARAMETERS - SEIMR/R-S (READ) 

SEI3RD 

Parameter 

SEIMR/R-S 

Parameter 

Description 

Measure 

Unit 

So, Si, S 2 , S 3 

Sst 

DELXag,st 

Probability of lo, Ii, h, h, ... recovering without worsening the clinical condition. 
Depends of the epidemic state and the age. 


TTo, Til, 712, 7^3 

Tlst 

PHIXag,st 

Time a patient in lo, Ii, h, I 3 , ... recovers 

day 

Tjo, Tji, r| 2 , r|3 

Tlst 

ETAXag,st 

Time a patient in lo, Ii, I 2 , 13 , ... to next infected state (develops symptoms, requires 
hospitalization, requires UCI, dies) 

day 


BIOLOGICAL PARAMETERS (CALCULATED) 

Para 

meter 

Equation 

SEIMR/R-S 

Parameter 

SEIMR/R-S 

Equation 

Description 

Measure 

Unit 

Yst 

l/7Ist 

GAMXag,st 

l/PHIXag,st 

Fraction of people in lo, Ii, I 2 , 13 ,... who recover 
in one day 

fpo/day 

ast 

l/r|st 

RHOXrg,ss 

l/ETAXag,st 

Fraction of people in lo, Ii, I 2 , I 3 , ... who 
complicate in one day 

fpo/day 


1/k 

FHII 

1/KAPP 

Inverse virus latency period 

1 /day 

AUXILIARY PARAMETERS 

S/st 

Sst 7st 

DEGXag,st 

DELXag,st X GAMXag,st 

S7st = Sst 7st 


crSst 

Sst (Jst 

RODXag,st 

DELXag,st X RHOXrg,ss 

ScTst = (1 -Sst) CTst = CTst - Sst CTst 


ScTst 

(1 -Sst) <Tst 

DERXag,st 

RHOXrg,ss - RODXag,st 

ScTst = (1 -Sst) CTst = CTst - Sst CTst = CTst " crSst 


SSst 

Sst 7st+ (1-Sst) CTst 

DGRXag,st 

DEGXag,st "1" DERXag,st 

S5st = (Sst 7st+ (1-Sst) CTst) = Syst 



The equations of SEI3RD model are the same as those of SEIR model considering the following changes: 

■ Equations 5S(t)/0t is modified 

5S(t)/at = - P X [ SsteINF Ist(t) ] X S(t) + X NP(t) - X S(t) 

■ Equations 8E{t)/dt is modified 
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5E(t)/5t = p X [ Sst.iNF Ist(t) ] X S(t) - \|/ X E(t) + k^x NP(t) 

■ Equation 0I(t)/5t is substituted by the equations 5Io(t)/5t 

01 o(t)/at = V X E(t) - 880 X Io(t) + A,' X NP(t) 

88st — (8st X yst + (l~8st) x c7st) 

■ Equation 0I(t)/5t is substituted by the equations 5Ist(t)/5t 

5Ist(t)/5t = 8c7st-i X Ist-i(t) - 88st X Ist(t) 

8ast = (1 -8st) X ast 

■ Equations 5R(t)/0t is modified 

5R(t)/5t = Zst.iNF 8yst x Ist(t) - |a'^ x R(t) + x NP(t) 


Syst — 8st X yst 


■ Equations 5D(t)/5t is modified 


5D(t)/5t = Sst.iF ^ X Ist(t) 
= 8on = (1 - 8n) on 


These equations are the equations of the SEIMR model 

3. SEIMR/R-S GENERAL EPIDEMIC MODEL 

The equations of the SEIMR model are presented in the next table. 



SEIMR - Differential Equations 

State 

St e SET 

State 

Increment 

State 

Decrement 

Natural 

Dead 

Exogenous 

Increment 

5S(t)/at 

su 

+ ^1= xNP(t) 

p X IS(t) X S(t) 

X S(t) 

X NP(t) 

5E(t)/5t 

EX 

P X IS(t) X S(t) 

¥ X E(t) 


X 

LU 

ei(t)/5t 

IN 

P X IS(t) X S(t) 



X' X NP(t) 

eio(t)/et 

10 

¥ X E(t) 

55st X Ist(t) 


X' X NP(t) 

5Ist(t)/5t 

IlF 

Sgsm X Ist-i(t) 




5R(t)/5t 

RE 

SstellF ^/st X Ist(t) 


n'' X R(t) 

X'^ X NP(t) 

5D(t)/5t 

ED 

Sst.IF c X Ist(t) 






AUXILIARY EQUATION 





IS(t) = [ SsteINF Ist(t) ] 




These equations serve to represent any of the three models studied. The conditions are as follows: 

1. SIR Model: 

■ Only consider one infected state IN the definitions of the infected state sets are: 

SteIF={IN}, steIlF={}, steI0={}, steINF={IN} and st£lN={IN} 

■ Do not include the exposed state (E), that means that steEX={} 

2. SEIR Model 

■ Only consider one infected state IN the definitions of the infected state sets are equal to SIR 
model. 

3. SEIMR Model 
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■ Considered multiples infected stats {10, II, 12,..., IN} IN associate to In, do not include the state 
lo that means: 

SteIF={lN}, SteIlF={Ii, I 2 , ... , In}, SteIO={Io}, SteINF={Io, Ii, I 2 , ... , In} and SteIN={} 

The next table presents the SETs of epidemic states needed to model the three epidemics models. The 
infected sets permit to model any of the three models with the same equations; they are used to define 
the existence conditions of the equations and in summation limits into the equations. 


EpicJemic States SETs 

Mo(jel 

Epidemic States 

Non Infected States 

Infected States 

STA 

SU 

EX 

RE 

ED 

ND 

INF 

IN 

10 

IlF 

IF 

SIR 

S, I, R, D, N 

S 


R 

D 

N 

I 

I 




SEIR 

S, E, I, R, D, N 

S 

E 

R 

D 

N 

I 

I 




SEIMR 

S, E, lo, Ii, I2, ... , In, D, N 

S 

E 

R 

D 

N 

lo, Ii, I2, ... , In 


lo 

Ii, I2, ... In-1 

In 


3.1. GENERAL FRAMEWORK 

To formulate the regional model segmented socio-demographically the following hypotheses are 
assumed: 

3.1.1. REGIONAL MODELING 

■ There is no contagion between people living in different regions. This can be true for large regions 
such as states or departments. But it is questionable for metropolitan areas (cities and conurbed 
regions) where there is intense traffic between regions. 

■ The inter-region interrelationship is modeled on the following assumptions: 

■ There is traffic of people between regions, which sets for each pair of regions the fraction of 
each segment, (pro,rg,ss, moving from the source region (ro) to the destination region (rg). 

■ In addition, the fraction of the time, (|)ro,rg,ss, is known to people from the region origin in the 
destination region during the period (one day). 

The calculation process involves determining the impact on the spread of the virus that the population 
movement has for this purpose it is calculated using the number of infected people who can move 
between two regions multiplied by the fraction of the time spent in the destination locality. This 
implies the following effects on the diffusion rate: 

1. Infected Movements 

■ Increasing the rate of diffusion in the destination locality due to those infected by coming from 
other regions, it is calculated as: 


Ilrg(t) — ZssESSR(rg) SroEROR(rg) (pro,rg,ss X (|)ro,rg,ss X ISro, ss (t) 

■ Decreased diffusion rate in the source region due to the infected by moving to other regions, it 
is calculated as: 


lErg(t) — ZssESSR(rg) ZrdERDE(rg) (prg,rd,ss X (|)rg,rd,ss X lSrg,ss(t) 


The net effect on the rg-region will be: 

lSrg,ss(t) = EsIeINF lst,rg,ss(t) 
IXrg(t) = ZssESSR(rg) lSrg,ss(t) 

IRrg(t) = IXrg(t) + Ilrg(t) ■ lErg(t) 


14 





















The SEIMR/R-S General Epidemic. 

Theory Validation and Application in Mexican Regions 




DECISIONWARE 


where 

Ist,rg,ss(t) fraction of the population infected in st-epidemic-state living in rg-region and ss- 


segment. 


ISrg,ss(t) fraction of the population infected living in rg-region and ss-segment. 

IXrg(t) fraction of the population infected living in rg-region 

Ilrg(t) weighted fraction of the population infected traveling to rg-region 

lErg(t) weighted fraction of the population infected traveling from rg-region 

2. Susceptible Movements 

■ Increasing the rate of diffusion in the destination locality (rg-region) due to those susceptible 
people by coming from other regions that may be infected in rg-region, it is calculated as: 


SIro,rg,ss (t) — (pro,rg,ss X (|>ro,rg,ss X Sro,ss(t) 


■ Decreased diffusion rate in the source region (rg-region) due to the susceptible people by moving 
to other regions that cannot be infected in rg-region, it is calculated as: 


SErg ,rd,ss(t) — (prg,rd,ss X (|)rg,rd,ss X Srd,ss(t) 


The susceptible population living in the rg-region ss-segment must be decremented by the 
susceptible people belonging to other regions: 


SNrg,ss(t) — Srg,ss(t) ~ Zrd€RDE(rg) SErg,rd,ss(t) 


where 

Srg,ss(t) fraction of the susceptible population living in rg-region and ss-segment. 

SIro,rg,ss(t) fractlon of the ss-segment susceptible population traveling from ro-region to rg- 
region 

SErg,rd,ss(t) fraction of the ss-segment susceptible population traveling from rg-region to rd- 
region 

3.1.2. SOCIO DEMOGRAPHIC SEGMENT MODELING 

■ An infected person in any segment can infect anyone susceptible in any socio-demographic 
segment. 

The calculation process implies that at the level of one region the population of any segment can 
infect the population of any other segment. To do this, diffusion (infection of susceptible from 
infected) is managed at a detailed level in the differential equations for all infected states, but in 
the differential equation the infected transmission is calculated based on the summation of all 
infected states in all ss-segments. 

The contagion of the susceptible population living in the rg-region and belonging to the ss-segment will 
be the sum of the contagions that occur in the rg-region (people that do not travel out of the rg-region) 
plus the contagions that occur in the rd-regions (destination regions) this is 


S2Irg,ss(t) — P X ( IRrg(t) X SNrg,ss(t) + Erd€RDE(rg) IRrd(t) X SErg,rd,ss(t) ) 


3.2. GENERAL FRAMEWORK 


The differential equations of the regional-segmented model are: 


5Srg,ss(t)/5t — - S2Irg,ss(t) - x Srg,ss(t) + ^^rg,ss x NP(t) 


5Erg,ss(t)/5t — S2Irg,ss(t) - \|/ X Erg,ss(t) -I" ^^rg,ss X NP(t) 
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St€lN 

St€lO 

StellF 

st€lF 


5Ist,rg,ss(t)/5t = S2Irg,ss(t) - yl^st,ss X Ist,rg,ss(t) + X^rg,ss X NP(t) 
5Ist,rg,ss(t)/5t = S2Irg,ss " yCJst,ss X Irg,ss(t) + X^rg,ss X NP(t) 
5Ist,rg,ss(t)/5t = yast-l,ss X Ist-l,rg,ss(t) - yast,ss x Ist,rg,ss(t) 
5Ist,rg,ss(t)/5t = yast-l,ss X Ist-l,rg,ss(t) - y(!^st,ss x Ist,rg,ss(t) 

5Rrg,ss(t)/5t = SstellF yCJst-l,ss X Ist,rg,ss(t) - X Rrg,ss(t) + SsseSSR(rg) X^rg,ss X NP(t) 


5Drg,ss (t)/5t — ZstellF |Liass X Ist,rg,ss(t) 
aNRrg(t)/at = X SRrg,ss(t) + x RRrg(t) 

Where y^s^ss, luass and yast,ss are rates (parameters) that depends on the configuration of the 
socio-demographics segments. In this case they are defined as 


y(^st,ss = I^ageAGS(ss) (yast,ag + |Llag) 

|Liass = I^ageAGS(ss) |Llag 
yCJst,ss = I^ageAGS(ss) yoCstag 

The definition equations of the regional-segmented model are: 

ISrg,ss(t) = ZsteINF Ist,rg,ss(t) 

IXrg(t) = SsseSSR(rg) ISrg,ss(t) 

Ilrg(t) — 2sseSSR(rg) 2roeROR(rg) (pro,rg,ss x (|)ro,rg,ss x ISro,ss(t) 
lErg(t) — I^sseSSR(rg) I^rdERDE(rg) cp''g/i'd,ss X (|)rg,rd,ss X ISrg,ss(t) 
IRrg(t) = IXrg(t) + Ilrg(t) - lErg(t) 

SRrg(t) = ZsseSSR(rg) Srg,ss(t) 

SIro,rg,ss(t) = (pro,rg,ss X (|)ro,rg,ss X Sro,ss(t) 
SErg,rd,ss(t) = (prg,rd,ss X (|)rg,rd,ss X Srd,ss(t) 
SNrg,ss(t) = S rg,ss(t) - Srd£RDE(rg) SErg,rd,ss(t) 
SINrg(t) = IRrg(t) x SNrg,ss(t) 

SIErg(t) — ZrdeRDE(rg) IRrd(t) X SErg,rd,ss(t) 
S2Irg,ss(t) — P X SINrg(t) + P X SIErg(t) 


RRrg(t) — SsseSSR(rg) Rrg,ss(t) 


DRrg(t) — EssESSR(rg) Drg,ss(t) 

From now on, the above mathematical definitions will be summarized as 

■C S, E, 1st, D, N } e © 

The next table shows the equations dividing the increment and the decrement on each state, it must be 
considered in the implementation of the mathematical models. The table includes the sets that defined 
the existence of the equations manly for the infected states. 



SIR Regional - Segmented Model - Differential Equations 

Set 

state 

state 

Increment 

state 

Decrement 

Natural 

Dead 

Exogenous 

Increment 


REGIONAL - SEGMENT EQUATIONS 

su 

5Srg,ss(t)/ dt 


S2Irg,ss(t) 

ILl'^ X Srg,ss(t) 

A.^rg,ss X NP(t) 

EX 

5Er9,ss(t)/5t 

S2Irg,ss(t) 

V|/ X Erg,ss(t) 


A.^rg,ss X NP(t) 

10 

5Ist,rg,ss(t)/ dt 

V|/ X Erg,ss(t) 

Y^st,ss X Ist,rg,ss(t) 


A.Vg,ss X NP(t) 
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SIR Regional - Segmented Model - Differential Equations 

Set 

State 

State 

Increment 

State 

Decrement 

Natural 

Dead 

Exogenous 

Increment 

IlF 


Y<^st-l,ss X Ist-l,rg,ss(t) 




IF 

Y<7st,ss X Ist,rg,ss(t) 



IN 

" S2Irg,ss(t) 


A-Vg,SS X NP(t) 

RE 

^Rrg,ss(t)/5t 

SlstellF YCJst,ss X Ist,rg,ss(t) 


X Rrg,ss(t) 

X NP(t) 

ED 

5Drg,ss(t)/ 5t 

|LXGss X Irg,ss(t) 




ND 

eNRr9(t)/5t 

X ( SRrg(t) + RRrg(t) ) 





SUSCEPTIBLE STATES EQUATIONS 


SRrg(t) = HsseSSR(rg) Srg,ss(t) 

SIro,rg,ss(t) = (pro,rg,ss X (t)ro,rg,ss X Sro,ss(t) 

SErg,rd,ss(t) = (prg,rd,ss X (|)rg,rd,ss X Srd,ss(t) 

SNrg,ss(t) = Srg,ss(t) “ SrdeRDE(rg) SErg,rd,ss(t) 

SINrg(t) = IRrg(t) X SNrg,ss(t) 

SIErg(t) = SrdeRDE(rg) IRrd(t) X SErg,rd,ss(t) 

S2Irg,ss(t) = 3 X SINrg(t) + 3 X SIErg(t) 

INFECTED STATES EQUATIONS 

ISrg,ss(t) = XstelNF Ist,rg,ss(t) 

IXrg(t) = IlsseSSR(rg) ISrg,ss(t) 

Ilrg(t) = SsseSSRCrg) 2roeROR(rg) (pro,rg,ss X (t)ro,rg,ss X ISro,ss(t) 

lErg(t) = SsseSSR(rg) SrdeRDE(rg) (prg,rd,ss X (t)rg,rd,ss X ISrg,ss(t) 

OTHER EPIDEMIC STATES EQUATIONS 

RRrg(t) = SsseSSR(rg) Rrg,ss(t) 

DRrg(t) = SsseSSRCrg) Drg,ss(t) 


3.3. DISCRETE EPIDEMIC MODEL 


The differential equations must be redefined in algebraic terms (discrete equations). The discrete 
equivalent equations of differential equations are 

[Srg,ss(t+At) - Srg,ss(t)]/At = - S2Irg,ss(t) - ILl*^ X Srg,ss(t) + X^rg,ss X NP(t) 

[Erg,ss(t+At) " Erg,ss(t)]/At = S2Irg,ss(t) - \|/ X Erg,ss(t) + A^rg,ss X NP(t) 

StelN 

[Ist,rg,ss(t+At) - Ist,rg,ss(t)]/At = S2Irg,ss(t) - yC^st,ss X Ist,rg,ss(t) + X^rg,ss X NP(t) 

StelO 

[Ist,rg,ss(t+At) - Ist,rg,ss(t)]/At = p X IRrg,ss(t) X Srg,ss(t) - yast,ss X Irg,ss(t) + X^rg,ss X NP(t) 

StellF 

[Ist,rg,ss(t+At) - Ist,rg,ss(t)]/At = yast-l,ss X Ist-l,rg,ss(t) - yast,ss X Ist,rg,ss(t) + X^rg,ss X NP(t) 

StelF 

[Ist,rg,ss(t+At) - Ist,rg,ss(t)]/At = yast-l,ss X Ist-l,rg,ss(t) - yl^st,ss X Ist,rg,ss(t) 

[Rrg,ss(t+At) - Rrg,ss(t)]/At = SstellF yCJst-l,ss X Ist,rg,ss(t) - ILX*^ X Rrg,ss(t) + SsseSSR(rg) X^rg,ss X NP(t) 

[Drg,ss (t+At) - Drg,ss(t)]/At = ZstellF |Llst,ag=AGS(ss) X Ist,rg,ss(t) 

[NRrg(t+At) - NRrg,ss(t)]/At = ILX*^ X SRrg,ss(t) + ILl*^ X RRrg(t) 


If At = 1 (one day), the above equations correspond directly to those that must be included in the 
optimization model using mathematical programming. 

Srg,ss(t+1) - Srg,ss(t) = - S2Irg,ss(t) - X Srg,ss(t) + X^rg,ss X NP(t) 

Erg,ss(t+1) - Erg,ss(t) = S2Irg,ss(t) - \j/ X Erg,ss(t) + X^rg,ss X NP(t) 
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StelN 


Ist,rg,ss(t+1) - Ist,rg,ss(t) — - S2Irg,ss(t) - yC^st,ss X Ist,rg,ss(t) + X^rg,ss X NP(t) 


StElO 


Ist,rg,ss(t+1) - Ist,rg,ss(t) — p X IRrg,ss(t) X Srg,ss(t) - yast,ss X Irg,ss(t) + X^rg,ss X NP(t) 


stEllF 



stElF 


Ist,rg,ss(t+1) - Ist,rg,ss(t) — yast-l,ss x Ist-l,rg,ss(t) - y(!^st,ss x Ist,rg,ss(t) 

Rrg,ss (t+i) - Rrg,ss(t) - SstellF yast,ag=AGS(ss) X Ist,rg,ss(t) - |Ll'^ X Rrg,ss(t) + ZsseSSR(rg) X^rg,ss X NP(t) 
Drg,ss (t+i) - Drg,ss(t) — SstsIlF |Jst,ag=AGS(ss) X Ist,rg,ss(t) 

NRrg(t+l) - NRrg,ss(t) = x SRrg,ss(t) + x RRrg(t) 

In addition, the definition restrictions should be considered 

■C S, E, 1st, D, N } e © 
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